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Abstract 

The multiple scatterings incurred by a hard quark produced in a nuclear 
medium induce the emission of soft gluons which carry a fraction of the quark 
energy and eventually affect the hadronization process. Here, the depletion of 
semi-inclusive hadron spectra in DIS on various nuclei (N, Ne, Cu, Kr) is com- 
puted as a function of v and z to leading order in a s through medium-modified 
fragmentation functions. Using the transport coefficient q previously determined 
from Drell-Yan production, the predictions are found to be in good agreement 
with EMC and HERMES preliminary data. Calculations on Xe targets are also 
presented and discussed. 



1 Introduction 

The interest of semi-inclusive hadron production in Deep Inelastic Scattering (DIS) on 
nuclear targets is at least twofold. First, and on rather general grounds, this reaction 
should help us to understand how a given QCD medium affects the non-perturbative 
hadronization mechanism. In particular, the use of nuclei — with a size R of a few 
fm — provides a direct information on the soft time scales A ~ 1/R involved in the 
dynamics of hadronization. In the second place, such a process also allows one to gain 
a quantitative insight into the scattering properties of cold nuclear matter. These may 
serve as a baseline to which properties of hot QCD media should be contrasted. This 
is clearly an important issue in high energy heavy ion physics where the quark gluon 
plasma is expected to be formed. 

Measurements of semi-inclusive hadron electroproduction on nuclei were performed 
a decade ago by the European Muon Collaboration (EMC) experiment who reported a 
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significant attenuation of charged hadron spectra in Copper as compared to Deuterium 
targets |T]. This was a strong indication that quark fragmentation functions get mod- 
ified in a nuclear environment 1 . More recently, the HERMES collaboration at DESY 
reported on high statistics measurements of charged hadron (h^), pion (tt^), kaon 
(K^, proton and antiproton (p, p) electroproduction on several nuclei (Deuterium, 
Nitrogen, Neon, Krypton) and on a wide kinematic acceptance jHlHj- These two data 
sets make quantitative analyses possible. 

Several suggestions have been advanced so far to account for the experimental 
results. Models based on the nuclear absorption of the produced hadron [oTIoTITllS] , gluon 
bremsstrahlung or partial deconfinement in nuclei [TUIIT]. were shown to describe 
the trend of the data fairly well. Following Wang and Wang ^T], we shall assume 
here that the observed hadron quenching actually comes from the energy loss incurred 
by hard quarks traveling through the nuclear medium. Although the effects of quark 
energy loss in DIS on nuclear targets have already been investigated in Refs. [TTlll2pi3j . 
a more complete analysis of EMC and HERMES data is still lacking. This is the aim 
of the present study. 

The outline of the paper is as follows. Hadron electroproduction in nuclei is com- 
puted to leading order in the coupling constant (Section^ using nuclear fragmentation 
functions modeled in Section El Predictions on hadron attenuation are then system- 
atically compared to EMC and HERMES preliminary data as a function of v and 
z (Section HJ). Calculations in Xenon targets in the HERMES kinematic acceptance 
will also be presented. Finally, we shall discuss and summarize our results in the last 
Section. 



2 Hadron production in DIS on nuclei 

Hadron production in semi-inclusive DIS can be computed in perturbation theory 
within the QCD improved parton model. To leading order (LO) in the strong cou- 
pling constant a s , the multiplicity N\ of a specific hadron species h in lepton- nucleus 
reactions is expressed in terms of nuclear parton densities f A (fq) and the q — > h 
(q — > h) fragmentation functions, D 1 ! (D^). Normalizing it to the number of DIS 
inclusive events N\ , it is given by 



1 dN\(y,z) 
N e A dv dz 



f dx^e 2 q [Zf^ A (x,Q 2 ) + (A-Z)f^ A (x,Q 2 )} a^(x,u) D h q (z,Q 2 , A) 
J dx^e] [Zf^ A (x,Q 2 ) + (A-Z)f-/ A (x,Q 2 )] a'^(x,u) (1) 



lr The attenuation of hadron electroproduction on a Carbon target was actually first observed at 
SLAC |2J. However, the SLAC measurements were not normalized to the number of inclusive events 
N% (see Eq. lfl]l) but the atomic mass number, hence rather reflected the nuclear modifications of 
parton densities (see discussion in Section |2J). For this reason, these data are not considered in the 
present study. 
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where v and Q are respectively the energy and the virtuality of the hard photon, and 
z = Eh/u is the energy fraction of the virtual photon carried away by the leading 
hadron. The nucleus is characterized by its atomic mass number A and its number of 
protons Z. The LO 7*g cross section appearing in (JTJ) is given by 
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where a em is the fine structure constant, M the nucleon mass, and s the center-of-mass 
energy squared of the lepton-nucleon reaction. The integral over Bjorken x = Q 2 jlMv 
in Eq. ((TJ is given by the Q 2 acceptance of the experiment. 

At small x < 0.05 and low Q 2 , the parton distribution functions (PDF) in nuclei 
fq differ significantly to those measured in a free proton, f? [T3]. These so-called 
shadowing corrections may therefore affect the hadron production cross section (JTJ. In 
the following, however, we shall not consider such corrections and assume fq^ A = f p 
for two reasons. First, the data we would like to discuss cover a Bjorken x range, 
x ~ 0.1, where the effects of nuclear shadowing prove rather small. Furthermore, the 
normalized hadron multiplicity Eq. JIJ is rather insensitive to the absolute magnitude 
of the parton densities. As noted in [7j, the shadowing corrections mostly cancel out 
in the ratio (JIJ) and finally affect hadron production by a percent at most. 

The parton densities in a proton f p we use in the present analysis are given by the 
MRST 2001 LO parameterization [T^j. It has been checked however that similar results 
were obtained using either the GRV LO ^H] and CTEQ5L [17] parton distribution 
functions. The neutron PDF are directly deduced from f p by isospin conjugation, 



i.e. vs 



d n , d p = u n , u p = d n , d p = u n , and s p 



In this paper, we would like to discuss the nuclear dependence of semi-inclusive 
hadron production. In particular, the theoretical calculations will be compared to the 
hadron production ratio in a heavy nucleus A with respect to a (light) Deuterium 
target, 



1 dN\{z,v) / 1 dN^(z,v) 



dv dz 



N 



D 



dv dz 



(3) 



which has been reported experimentally. When x is not too small, hadron production 
(0) is dominated by the fragmentation of valence quarks. Assuming for simplicity that 
only the up quark channel contributes to the process, the double production ratio 
Eq. (JHJ) approximately reduces to 



R h A (z,u)^D h u (z,Q 2 ,A) / D h u (z,Q 2 ,D). 



(4) 



Therefore, hadron production in DIS on nuclei is mostly sensitive to the nuclear depen- 
dence of the fragmentation functions. In particular, the reported depletion of hadron 
multiplicity on heavy nuclei R\(z, v) < 1 on a wide kinematic range is a clear hint that 
fragmentation functions are somewhat modified in the medium. 
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The starting point to determine the hadron attenuation R\{z, v) is therefore the 
computation of nuclear fragmentation functions D(z,Q 2 ,A). This task is carried out 
in the following section. 



3 Nuclear fragmentation functions 

3.1 Model for medium- modified fragmentation functions 

To leading order in the coupling and at leading twist, the virtual photon picks up a 
quark 2 in the nucleus which subsequently hadronizes into the observed leading hadron. 
In presence of a QCD medium, however, the hard quark suffers multiple scattering and 
radiates soft gluons all along its path. Due to this medium-induced gluon radiation, 
the quark energy is reduced from E = v to E = v — e at the time of hadronization. 
This picture is schematically depicted in Figure [TJ 



A 




Figure 1: Schematic representation of the semi-inclusive hadron electroproduction in a 
nuclear target. The virtual photon couples to a hard quark which subsequently emits 
soft gluons while escaping the nucleus. 

How to relate these final state interactions to nuclear fragmentation functions re- 
mains however unclear. Such an attempt has been performed recently in a series of 
papers [TTIITS] in which a higher-twist perturbative framework has been applied suc- 
cessfully to describe hadron production in DIS and heavy ion reactions. 

Here, we shall rather adopt the effective model suggested in Ref. Within this 
approach, the quark energy shift leads to a rescaling of the momentum fraction in 

2 In the following discussion, we shall ignore the fragmentation of antiquaries as x is not too small. 
The theoretical computations Eq. consider of course this channel as well. 
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presence of a QCD medium, 

Eh * E h 



- T-t-T, (5) 
v v — e 1 — e/v 

where Eh is the measured hadron energy and e the energy lost by the hard quark 
while going through the nucleus. Consequently, the fragmentation functions in nu- 
clei Dq(z,Q 2 , A) may simply be expressed as a function of the standard (vacuum) 
fragmentation functions D q (z,Q 2 ) through [TU] 

rv-E h 

z D h q {z, Q 2 , A) — de D(e, v) z* D h q (z\ Q 2 ). (6) 

Jo 

Here, D(e, v) denotes the probability for a quark with energy E — v to lose an energy 
e. This probability distribution, or quenching weight, will be made explicit hereafter. 
Since fragmentation functions fall steeply at large z, even the small shift Az = z* —z ~ 
zejv in JHJ) may substantially suppress the hadron yields in nuclear targets. 

Note that in the Bjorken limit v ^> e, hence Q ^> e at finite x, the effects of the 
final state interactions become negligible, z ~ z* . No more nuclear dependence of 
the fragmentation functions is observed D q (z,Q 2 , A) = D q (z,Q 2 ), in agreement with 
QCD factorization theorems 



The vacuum fragmentation functions D q (z, Q 2 ) into various hadron species (charged 
hadrons, pions, kaons) which enter Eq. (JBJ) are here given by the leading order analysis 
of e + e~ data performed by Kretzer [2T] . 

Let us now turn to the computation of the quenching weights for hard quarks going 
through nuclear matter. 



3.2 Quenching weight 



As apparent in Eq. (JBJ), the knowledge of the probability distribution D(e, v) is essential 
to determine the nuclear fragmentation functions. This point has been raised recently 
by Baier, Dokshitzer, Mueller, and Schiff (BDMS) in Ref. [22] i n which they express 
Die, u) as a Poisson series 



n=0 



i=l 



/ dwi 



dJjujj, v) 



i=l 



(7) 



assuming gluon emissions to be independent. 



The expression Eq. J7J) explicitly relates the probability distribution in the energy 
loss to the gluon spectrum dl /du radiated by hard quarks produced in QCD media, 
which has been computed perturbatively by Baier, Dokshitzer, Mueller, Peigne, and 
Schiff (BDMPS). To first order in the quark energy O it reads j2H] 



where Cp = 4/3 is the Casimir operator in the fundamental representation and a s = 
g 2 /47r ~ 1/2 the strong coupling constant. Note that a s in (jSJ) is evaluated at a 
soft scale given by the virtuality of the radiated gluons, while the renormalization 
scale appearing in the parton densities and fragmentations functions in is given 
by the photon virtuality, Q. The appearance of these two scales follows our implicit 
assumption on the factorization between the production process and the final state 
interactions. The typical gluon energy in the spectrum is set by the relevant scale 

u c = \qL 2 (9) 

which characterizes the QCD medium. The gluon transport coefficient q reflects the 
medium gluon density (say, its scattering property) while L is the length of matter 
covered by the hard quark. Note that q can be very large in practice in a hot pion gas 
or quark gluon plasma [231 • Energy loss of partons may thus be used to probe such 
dense media expected to be formed in high energy heavy ion reactions (23H2H]- The 
absolute values for uj c in nuclei will be discussed below. 

Using Mellin transform techniques, BDMS resumed the Poisson series to eventually 
give the quenching weight Eq. (J7J) a simple integral representation. This allows for 
the numerical computation of -D(e, v) from the BDMPS gluon spectrum (jSJ) performed 
in [T2]. The distribution has been given a simple analytic parameterization which we 
shall use in the present analysis. 

Transport coefficient for nuclear matter q 

The quenching of hadron spectra depends crucially on the transport coefficient for 
nuclear matter. The higher the q, the larger the typical quark energy loss e oc q, 
hence the stronger the depletion in large nuclei. Based on the (not too) small x gluon 
distribution in a nucleus, BDMPS provide a perturbative estimate for the transport 
coefficient in nuclear matter, q ~ 0.25 GeV/fm 2 |26j . 

A slightly larger — although consistent — value has been obtained recently from 
our LO analysis of Drell-Yan production data [2Z] measured in n~-A collisions by the 
NA3 collaboration, 

q = 0.72 GeV/fm 2 . (10) 

This value will be taken for the numerical applications to come. It should already 
be kept in mind, however, that this value could not be determined with a high ac- 
curacy (la = 0.54 GeV/fm 2 ) because of the large statistical errors in the Drell-Yan 
measurements. 

As we shall discuss later (Section [S]), hadron depletion in DIS on light nuclei (say, 
Nitrogen and Neon targets) may provide an independent estimate for the transport co- 
efficient q. In particular, the high statistics HERMES data should prove of considerable 
help to constrain further the estimate (|10j) from Drell-Yan production. 

Length L 

Once q is fixed, only the mean length L covered by the hard quark in the nuclear 
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medium remains to be determined. Assuming that hadronization occurs well after the 
hard quark has escaped from the nucleus, in other words when the hadron formation 
time tf is large, L may simply be estimated for any given density profile p(r). Taking 
for simplicity a sharp sphere density p(r) = p x 0(f?, — |r|), it is given by 

3 

L max = - R, (11) 

where the nuclear radius R = (3/(47rp)) 1//3 A 1 / 3 ~ 1.17 A 1 / 3 fm with p = 0.15 fm~ 3 . 

When tf gets smaller than twice the nuclear radius, however, quarks produced on 
the "back" side of the nucleus hadronize while still inside the medium. The length (jllj) 
given above is therefore no longer correct. Integrating over the nuclear volume, we now 
obtain 



L(t f ) =t f x 



3*/ , 1 /V 



8 R 64 V R 



if t f < 2 R. (12) 



Of course, L(tf) = L max at tf = 2R. For small formation times (large nuclei), tf/R<^ 
1, energy loss effects are strongly reduced as the quark has propagated over shorter 
distances than (fTTj) . L(tf <C R) — tf <C L max . It is implicitly assumed here that 
a color singlet hadron (Ca = in (jHJ)) does not lose any energy by gluon radiation. 
The produced hadron in the medium may nevertheless interact and be absorbed in 
the nucleus. In this approach — that we want to keep as simple as possible — , this 
mechanism is not taken into account. However, we shall come back to this point in the 
discussion section. 

The soft scales involved in the hadronization process prevent any perturbative cal- 
culation of the formation time tf. Therefore, the absolute value for tf as well as its 
kinematic dependence is not clearly known. We shall take the Bialas-Gyulassy esti- 
mate [B] obtained in the framework of the Lund model 3 

'/~( mi/ ^-; + *> t- < i3 > 

where a is a non-perturbative scale in the order of the string tension, although slightly 
reduced in a nuclear environment [7]. We take a = 0.75 GeV/fm throughout this 
study. The formation time (|T3~|) exhibits a maximum around z ~ 0.3 and decreases like 
tf ~ (1 — z) at large z, observed as well within gluon bremsstrahlung models 0. 

Plugging the Lund formation time (|T3|) into Eq. ()12)) . the mean (normalized) length 
L/R is plotted as a function of the quark energy v and the hadron momentum fraction 
z in Figure El Also shown as a dashed line is the asymptotic length L max / R when 
hadronization occurs outside of the nucleus. We observe that L gets much smaller 
than L max when z approaches or 1, as tf vanishes. We shall thus expect in these 
specific kinematic limits a dramatic discrepancy with the prediction L = L max assuming 
hadronization outside of the nucleus. Because of the Lorentz boost L oc v in (J13|) . 

3 It has been stressed in [7] that Eq. (|13|) slightly overestimates the Lund model predictions with 
standard parameters at small z. We checked however that this does not strongly affect our results. 
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Figure 2: Mean length (|12j) covered by the hard quark in the nucleus of radius R as 
a function of z and for different values of the quark energy v. The asymptotic length 
L max assuming hadronization outside of the nucleus is shown as a dashed line. 



highly energetic quarks v ^> oR will fragment into hadrons after escaping the nuclear 
medium. As discussed in the next section and already apparent in Figure |2l finite 
formation time effects become irrelevant to describe high energy data. 

To give the reader a feeling for such formation time effects, hadron attenuation will 
be computed in the following account assuming both the asymptotic length L max (fTTj) 
(dashed lines) as well as the more realistic estimate Eq. (jl2j) (solid lines). 



4 Results 



Hadron production in semi-inclusive DIS on nuclear targets may now be determined 
from the perturbative expression using the nuclear fragmentation functions de- 
rived in the latter section. Calculations are performed using the acceptance of the 
experiments. We have used the mean kinematic variables (u), (z), and (Q 2 ) measured 
experimentally whenever they were available, as it was shown in [7] that these were 
slightly underestimated (by less than 10%) by the LO prediction (JTJ). 

The quenching of hadron multiplicity R^iy) and R^z) is compared to the EMC £Q 
and HERMES [3] data on Cu as well as on N and Kr nuclei, respectively. We shall also 
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compare in this section our results to the HERMES preliminary data on Neon 0j. Fi- 
nally, predictions on a Xenon target in the HERMES kinematic window are presented. 



4.1 v dependence 



We start the section by first considering the quark energy dependence v of the hadron 
quenching. As previously emphasized, the energy loss effects should vanish as the 
energy v gets large as compared to the typical scale uj c of the medium. Naturally, the 
quenching factor R\(v) should therefore increase with v and eventually approaches 
unity at asymptotic energies. 

This is what is observed experimentally. In Figure H3 is plotted the charged hadron 
attenuation R h± (u) in a Cu nucleus on a wide range, 15 < v < 200 GeV. The calcula- 
tions are found to be in good agreement with the measurements performed by the EMC 
collaboration. Notice that these high energy data do not exhibit any formation time 
effects in such a small nucleus, vjaR > 6, as shown from the tiny difference between 
the dashed and solid lines in Figure 01 
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Figure 3: Quenching of charged hadron spectra in Cu target as a function of v. The 
theoretical calculations with (solid) and without (dashed) finite formation time effect 
are compared to EMC data pQ. 



Let us turn to lower energy data. The HERMES collaboration measured semi- 
inclusive hadron production on Nitrogen and Krypton targets. What is more, the 
experimental apparatus allows for a particle identification of various hadron species, 
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such as pions, kaons, protons and anti-protons. Since no fragmentation functions are 
provided in the (anti-)proton channel, only the quenching of charged hadron, pion, and 
kaon spectra will be presented. 

Figure |U (upper panel) indicates that the quenching of charged hadrons on Nitro- 
gen is extremely well reproduced by our predictions, while the calculations slightly 
overestimate the attenuation in the pion (isospin averaged) sector. It should be noted 
that the stronger pion quenching K*{y) < R h (v) reflects only the different acceptance 4 
in the HERMES experiment. In the Kr target, the suppression of charged hadrons is 
again well described by the calculations. The trend predicted in the pion channel in 
Kr is close to what is observed experimentally, although the absolute magnitude of the 
pion quenching slightly underestimates the data (~ 5 — 8%). 

Perhaps more interesting are the HERMES preliminary measurements on kaon 
quenching in krypton, (is) (figure El lower panel). A significantly larger deple- 
tion is observed in K~ production on the whole energy range, 
Furthermore, this noticeable isospin dependence in the kaon sector is apparent in the 
theoretical predictions as well. It is indeed clear from Eqs. (jlj) to (jUJ) that the quenching 
is directly related to the slope of the vacuum fragmentation functions 5 



where rj^ is the slope of the fragmentation functions at large z, D^(z, Q 2 ) ~ (1 — z) v ™. 
Schematically, the up quark will preferably fragment into a K + by picking up easily 
a surrounding s quark in the sea (valence- type fragmentation channel). Inversely, the 
difficulty to produce a negative kaon from a hard up quark (sea-type) will translate 
into a steeper slope, rj^ = + 1 and thus in a stronger quenching Eq. ((HJ). At very 
small a; < 1, semi-inclusive production originates from the sea quark fragmentation, 
u = u 6 . We therefore predict that no more isospin effect should be observed in this 
region, and thus R K (u) ~ R K (y) 

Formation time effects start to become visible in Krypton targets when v is not too 
large, v < 12 GeV (see Figure HJ). Nevertheless, the rather small differences between 
the dashed and solid curves indicate that the krypton nucleus (A=84) is not too large 
to probe the hadronization in the medium. More pronounced effects are expected in 
heavier nuclei. 

Finally, hadron attenuation has been computed on a Neon target (Figure upper 
panel) measured very recently by the HERMES collaboration While the observed 

4 More precisely, it comes from the different cuts in the z acceptance. Pions were measured with 
z > 0.5 while a cut z > 0.2 has been applied to the charged hadron channel. 
5 We omit, for clarity, the integral J deD(e) in Eq. I|14fl . 

furthermore, the sea is almost flavor symmetric u ~ d at small x. However, down quarks will be 
disfavored due to their electric charge. 



dD h u 1 z e 
dz D h u (z)~ 



1-VuX 



z e 



r/(l -z) 



(14) 
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Figure 4: Quenching of charged hadron, pion, and kaon spectra in N (upper) and 
Kr targets (lower) as a function of v. The theoretical calculations with (solid) and 
without (dashed) finite formation time effect are compared to HERMES data taken 
from Ref. The 4% systematic errors are not shown. 
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Figure 5: Quenching of charged hadron, pion, and kaon spectra in Ne (upper) and 
Xe targets (lower), with (solid) and without (dashed) finite formation time effect, as a 
function of v. The HERMES preliminary data on Ne nuclei are taken from Ref. [I]. 
The 3.3% systematic errors are not shown. 
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pion quenching is well reproduced by the calculations, the predicted attenuation notice- 
ably overestimates the preliminary data in the kaon sector, and particularly in the K~ 
channel. This feature was already present, although less pronounced, in the compari- 
son with Krypton measurements. The underestimate of the slopes of the flavor-tagged 
u — > K + (K~) fragmentation functions, poorly constrained by e + e~ data, may be at 
the origin of this discrepancy (23- m particular, the assumption made in Ref. |2I| that 
the sea type channels (e.g. D£ ) get suppressed by one extra power in (1 — z) (see 
discussion above) seems to be disfavored by experimental data on kaon production |28j . 

Figure El (lower panel) also displays our predictions on Xenon nuclei (A=131), in 
which formation time effects get larger. As later stressed in the discussion section, 
measurements in such a heavy nucleus will be of crucial importance to disentangle the 
possible various mechanisms responsible for the observed hadron attenuation. 

4.2 z dependence 

Another way to probe the dynamics of hadronization in a nuclear medium is to look 
at the z dependence of hadron quenching. In our calculations, we expect a stronger 
hadron quenching as z increases for two reasons. First, as can been seen from the rough 
estimate Eq. (|14j). the slope of the fragmentation functions increases with z, and thus 
so does the attenuation R\(z). 

The second reason that leads to significant hadron depletion at large z actually 
comes from the restricted phase space. Measuring hadrons with large momentum frac- 
tion z selects on those (rare) fragmenting quarks that have lost a tiny amount of en- 
ergy. Assuming for simplicity that the fragmentation functions are flat, z*D(z*, Q 2 ) ~ 
zD(z,Q 2 ), the quenching factor may be approximated by 



which vanishes as z approaches 1. Let us note however that the ratio (jl5|) will tend 
to a constant po > in the parton energy loss frameworks developed by Gyulassy, 
Levai, Vitev and Wiedemann |3U] i n which the quenching weight D(e, v) contains 
a discrete part, poS(e) Similar behavior is also observed within the BDMPS 

formalism when ultra soft gluons are removed from the perturbative spectrum |12j . 

Looking at Figure |Hl our predictions on the z dependence of the charged hadron 
quenching R h± (z) in Copper targets describe well the EMC measurements, with the 
exception of the highest z bin. Again, high energy quarks hadronize well after escaping 
the nucleus in this whole z range. 

On the contrary, we do observe huge formation time effects in the HERMES data 
which are taken at much lower energy (u) ~ 10-20 GeV (Figure |7j). The predictions 
fairly reproduce the HERMES measurements for all hadron species when z is not too 
high. At large z > 0.7, nevertheless, the discrepancy between the predictions assuming 
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Figure 6: Quenching of charged hadron spectra in Cu target as a function of z. The 
theoretical calculations with (solid) and without (dashed) finite formation time effect 
are compared to EMC data pQ. The dashed and solid lines are superimposed. 



L = I/max (dashed lines) and the data is striking. As discussed above, the hadron 
attenuation R(z) falls too steeply as z goes to one, unlike the experimental results. 
This does not come as a surprise as we have already anticipated in the latter section 
that this assumption is no longer justified since L ~ tf ~ (1 — z)vja <C L max . 

When the realistic length Eq. (fT2j) is used, the typical scale u c oc L 2 ~ (1 — z) 2 
is getting smaller and weakens the quark energy loss — and subsequently the hadron 
attenuation — as 2 increases. This effect somehow compensates the quenching due to 
the slope of the fragmentation functions (fT4"|) and to phase space restriction (fT3j) at 
moderate z, and eventually takes over at z > 0.9 (see Figure |7j). This non-monotonic 
behavior of the hadron depletion is actually a manifestation of the quadratic behavior 
of the parton energy loss in QCD media. We may indeed convince ourselves of this 
by replacing the energy loss e oc uj c ~ (1 — z) 2 in the rough estimate (J14)) . It is 
however doubtful whether this could seen in the data as the final state interaction of 
the produced hadron may somehow modify this picture. As shown in Figure [7J a nice 
agreement is reached between these calculations (solid lines) and the HERMES data 
for all z. 

The predicted hadron attenuation in Ne and Xe targets is also plotted in Figure |H1 
as a function of z. The preliminary HERMES measurements in Neon are well described 
by these predictions. Once more, let us emphasize that measurements in Xenon targets 
will prove essential to further constrain the hadron formation time estimates at large z. 
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Figure 7: Quenching of charged hadron, pion, and kaon spectra in N (upper) and 
Kr targets (lower) as a function of z. The theoretical calculations with (solid) and 
without (dashed) finite formation time effect are compared to HERMES data taken 
from Ref. The 4% systematic errors are not shown. 
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Figure 8: Quenching of charged hadron, pion, and kaon spectra in Ne (upper) and 
Xe targets (lower), with (solid) and without (dashed) finite formation time effect, as a 
function of z. The HERMES preliminary data on Ne nuclei are taken from Ref. [3]. 
The 3.3% systematic errors are not shown. 
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Let us finally stress that theoretical uncertainties might get large in the high z — > 
1 region. In particular, important higher twist corrections may come into play in 
this restricted region of phase space f2EI; given moreover the fact that the typical 
values Q 2 ~ 2 — 4 GeV 2 measured in these reactions are pretty small. Moreover, the 
resummation of soft and collinear gluons leads to large contributions a s log(l— z) / (1— z) 
in the fragmentation functions jH2] which are not taken into account in the leading order 
analysis we use here. 



5 Discussion and summary 

The nuclear dependence of semi-inclusive hadron production in DIS comes in our cal- 
culations from the sole effect of the quark energy loss in nuclei. Of course, many other 
mechanisms, which have not been considered here, could play a role in the hadron 
quenching observed experimentally. The absorption of the (pre)hadrons produced in 
the medium may for instance substantially affect its production 0ini0IE] • This will be 
particularly true at large z and/or small v when the formation time is small. It was 
also stressed that the partial deconfinement which could occur in large nuclei rescales 
the fragmentation functions and eventually suppress the hadron yields [T(Hl7]. 

However, although such effects have not been taken into account here, a very nice 
description of the EMC and more recent HERMES preliminary data for several hadron 
species and various nuclei is still achieved. The already good agreement assuming 
only energy loss effects in our approach is therefore a hint that the above mentioned 
competing processes play a rather modest role. Let us moderate this statement. The 
strength of each of these mechanisms is governed by one parameter which is in general 
poorly constrained. Take the quark energy loss for instance. We have used here the 
transport coefficient for nuclear matter q = 0.72 GeV/fm 2 extracted from Drell-Yan 
data, whose absolute magnitude is not known with a great accuracy. Similarly, the 
nuclear absorption strongly depends on the pre-hadron cross section a* about which 
little is known. 

It is therefore a delicate matter to disentangle these various processes. We suggest 
a couple of observables which may help us in this prospect. First, we have seen that the 
mean length L ~ tf, hence the scale uj c , becomes independent of R when the nuclear 
radius is large, R ^> vjo. The effects of quark energy loss on hadron quenching thus 
somehow saturates when the nuclei are sufficiently heavy. This can already be seen from 
the comparison of the rather similar hadron quenching in Kr and Xe targets, i?Kr — Rxe 
(Figures [7| and |Hl solid lines). It is clearly at variance with the absorption models 
which predict of course a much stronger attenuation as the nucleus gets larger. We 
therefore predict that the ratio of hadron depletion in two large nuclei (say, Xe/Kr) will 
provide a good measurement of hadron absorption, only. Inversely, hadron quenching 
in small nuclei could help to constrain the transport coefficient in nuclear matter, 
without being spoiled by hadronic absorption 7 . Fitting q to the charged hadron and 

Let us however mention that rescaling models predict a non-negligible suppression in such small 
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pion quenching in Nitrogen and Neon data, we found q = 0.59 ± 0.18 GeV/fm 2 and 
q = 0.70 ± 0.06 GeV/fm 2 , respectively, in excellent agreement with the Drell-Yan 
estimate 8 , Eq. (fTU|) . which suffered from a larger statistical uncertainty. 

We also think that the isospin dependence of the kaon quenching could be im- 
portant. As already mentioned, a similar attenuation ratio R K (is,z) ~ R K+ (v,z) is 
expected at small x if quark energy proves the most relevant mechanism at work. On 
the contrary, a stronger K~ quenching should still be observed in the framework of 
absorption models from the larger interaction of negative kaons. The Q 2 dependence 
(at a given v) of the ratio R K+ / R K may thus clarify the relevant underlying processes. 

To conclude, a leading order calculation of semi-inclusive hadron production in DIS 
on nuclei has been performed. The fragmentation functions were modified to take 
into account the energy loss incurred by hard quarks while they propagate through the 
nuclear medium. The quenching of hadron spectra was then compared to the EMC and 
HERMES preliminary data in several nuclear targets, as a function v and z. Using the 
transport coefficient determined from Drell-Yan production, a nice agreement between 
our predictions and the experimental measurements is observed for the various hadron 
species. In particular, we emphasized that hadron production in small nuclei may help 
us to determine and constrain what the transport coefficient (hence the quark energy 
loss) of nuclear matter is. Measurements of semi-inclusive hadron production in heavy 
nuclei are highly desirable since, in addition to be a sensitive probe of formation time 
effects, these prove crucial to disentangle quark energy loss on the one hand from 
hadron absorption on the other hand. 

Finally, it should be interesting to constrain, from the HERMES and EMC data, 
the rapidity dependence of the quenching of charged hadron and pion hadroproduction 
inclusive spectra in deuteron-gold collisions at RHIC energy [55] , 
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